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Abstract 

We discuss bosonization in three dimensions of an SU(N) massive 
Thirring model in the low-energy regime. We find that the bosonized 
theory is related (but not equal) to SU(N) Yang-Mills-Chern-Simons 
gauge theory. For free massive fermions bosonization leads, at low 
energies, to the pure SU(N) (level k = 1) Chern-Simons theory. 
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1 Introduction 



In this paper we investigate the problem of identifying a bosonic equivalent 
of a theory of self-interacting fermions with SU (N) symmetry in 2 + 1 di- 
mensions. In a previous publication |L], two of us showed that the low-energy 
sector of the U(l) massive Thirring model in 2 + 1 dimensions is equivalent 
to the Maxwell- Chern- Simons gauge theory. Here we extend that analysis to 
a theory with non-abelian SU(N) symmetry. We show that, just as in the 
abelian case, it is possible to bosonize the low-energy regime of the theory. 
It is also a gauge theory and, as expected, it is non-abelian and it is closely 
related to a (level k — 1) SU(N) Chern-Simons gauge theory. Naively, one 
would have expected that in the non-abelian case the bosonic theory should 
be the Yang-Mills-Chern-Simons (YMCS) theory. Surprisingly, we find that 
the bosonized theory is related, but not identical, to YMCS. However, in the 
limit of weak (fermionic-) coupling g — > 0, the bosonic theory becomes equal 
to the CS gauge theory. At first sight it may seem surprising that a simple 
theory such as free massive fermions may be equivalent to YMCS. It may be 
noted that also in 1 + 1 dimensions the bosonized theory of a massive free 
fermion is also non-trivial: it is equivalent to a special case of the sine-gordon 
theory P, |3|] (in the abelian case) and to a perturbed Wess-Zumino-Witten 
theory |§. 

We will follow the same strategy as in reference and seek a bosonic the- 
ory which reproduces correctly the low- energy regime of the massive fermionic 
theory. It should be stressed that this procedure is, in a sense, opposite to 
what is done in 1 + 1 dimensions. There, bosonization is a set of operator 
identities valid at length scales short compared with the Compton wavelength 
of the fermions. Here, instead, we only consider the long distance regime. 

This paper is organized as follows. In Section |2] we present the mapping of 
the low energy sector of the non-abelian fermionic theory into an equivalent 
gauge theory. In Section [3] we derive a set of identities for the fermionic 
currents. In Section [| we discuss the role of the Wilson loops of the gauge 
theory in the equivalent Fermi theory. In Section^] we draw a few conclusions 
on the mapping presnted here. 
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2 The Mapping 

We start from the three-dimensional (Euclidean) massive SU(N) Thirring 
model Lagrangian: 

C Th = j>(i$ + mW-^j a % (1) 

where ip is a two-component Dirac spinor in the fundamental representation 
of SU(N), and j a ^ the SU(N) current, 

3 ail = tf*$7V. (2) 

The coupling constant g 2 has dimensions of inverse mass. (Although non- 
renormalizable by power counting, four fermion interaction models in 2 + 1 
dimensions are known to be renormalizable in the 1/N expansion ||.) We 
normalize the SU(N) generators according to tr t a t b = 5 ab . 
The partition function for the theory is defined as 

2 

Z Th = I ZtyXty exp{- J d 3 x m$ + m)V> - yj^]}- (3) 

We eliminate the quartic interaction by introducing a vector field taking 
values in the Lie algebra of SU(N), through the identity 

exp{|- J d 3 x j a X) = J Va„ exp{- tr J d 3 x (^a% + fa^} (4) 

(up to a multiplicative normalization constant), so that the partition function 
becomes 

Z Th = J V^VijVa^exp{- J d 3 x + m + 4)^j + ^tra%]}. (5) 

The fermionic path-integral can now be done and it yields, as usual, the 
determinant of the Dirac operator 

J T>4>Vipexp{-Jd 3 xijj(i$ + m + <jl)?p} = det(i$ + m + $) (6) 

Being the Dirac operator unbounded, its determinant requires regulariza- 
tion. Any sensible regularization approach (for example, ^-function or Pauli- 
Villars) gives a parity violating contribution IJ-pf. There are also parity 
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conserving terms which have been computed as an expansion in inverse pow- 
ers of the fermion mass: 

lndet(z$ + m + <f) = ±-Ls CS [a] + Ipc[a] + 0(<9 2 /m 2 ), (7) 

107T 

where the Chern-Simons action Scs is given by 

Scs[a] = J d 3 xie^ x tr J d 3 x (f^ax - ^a^a u a x ). (8) 
Concerning the parity conserving contributions, one has 

Ipc[a] = -^—trf<& rf, v + ■ • • , (9) 

where 

fp> = 9^a v ~ 9 v an + [a M , a v \. (10) 
Using this result we can write Z^ h in the form 

Zth = J X>a jU exp(-5 e/ /[a]), (11) 

where S e ff is given by 

SefM = ^2 tr jAa/T^ cs H + ^ti j d 3 x T f, u + 0(d 2 /m 2 ). 

g vr vrm ^ 

Up to corrections of order 1/m, we recognize in S e jf the non-abelian version 
of the self-dual action Ssd introduced some time ago by Townsend, Pilch 
and van Nieuwenhuizen ||, 

Ssd [a] = [d 3 x^ tr T -^S CS [a]. (13) 

J Ig* 107T 

Then, to leading order in 1/m we have established the following identification: 

Z Th ~Z SD = J Va^ex P (-S SD [a]) (14) 

In the abelian case, it has been shown by Deser and Jackiw [HJ that the 
model with dynamics defined by Ssd is equivalent to the Maxwell- Chern- 
Simons (MCS) theory. Using this connection, we have shown in M the 
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equivalence of the abelian Thirring model and the MCS theory. Our proof 
was based in the use of an "interpolating Action" Si connecting the self 
dual and the MCS actions. It has been already recognized in [ID] that the 
non- abelian extension of these kind of equivalences is more involved. In fact, 
we shall show here that the non-abelian self-dual action (and, consequently, 
the SU (N) Thirring model) is not equivalent to a Yang-Mills-Chern-Simons 
theory (the natural extension of the abelian MCS theory) but to a model 
where, instead of the Yang-Mills action, one has a more complicated term, 
which reduces to the quadratic tr F 2 U only in the g 2 — > limit. 

To see this, we consider, following ||13||, the interpolating Lagrangian 



1 i 1 

L z [a,A] = ^tra^ ± — e^ A tr a^F^ + A u a x ) ± —-L C s[A], (15) 
Sg z p lo/T lo/T 

with F^v the non-Abelian curvature 

F^ = d^Av - d v A^ + [A^, A v \. (16) 

and L^s [A] is the Chern-Simons Lagrangian. The partition function associ- 
ated to Li is defined as 

Zi = J VA^Va^M-J d 3 xLi[a,A]}. (17) 

In order to see the connection between Zi and the partition function Zsd 
for the self-dual system, we shift in ( |PTD the A^ variables as follows: 

Afi — Afj, — dfj,- (18) 

Under this change, one has 

F„ V [A] = F^A] - D[A] [lx a u] + [a M , a v ], (19) 
with the covariant derivative defined as 

D ll [A]=d ll +[A (t , }. (20) 
Concerning the Chern-Simons Lagrangian, one has 

L CS [A] = L CS [A\ - 2e fMuX a fl (F uX [A] - D u [A]a x + ^a u a x ). (21) 
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Puting all this together we easily find 

2g 2 16tt GM J 8vr m ^ 3 



Li[a,A] = —tra^ ± -^-L CS {A] =f i— tre^ A a M (9„a A + -a u a x ). (22) 



Transformation (|18"D has unit Jacobian and completely decouples the A field. 
Denoting by N the result of integrating over A, one then has 

Zj = N J Pa M exp{-^tr Jd 3 xa^} x 

lr 2 
exp{±i — tr J d 3 xe^ uX a^(d u a x + -a u a x )}. (23) 

We recognize in this expression the action for the self-dual system defined in 
eqs. (|13"D. Then, we have proven that 

Zj = N j Va,exp(~S SD )- (24) 



Comparing eqs. (H|) and ( |24|) we can establish the following relation 

Z Th « Zj. (25) 

We shall now proceed to perform the path-integrations in Zj in the inverse 
order, that is, first integrating over a M . To this end, starting again from 
eq. fll5| ) we write 

Z I = J VA^expiT^ScslA}) x exp(-/[A]), (26) 

with 

exp(— I[A]) = J Va fl exp{ — -^tr JcPxa^} x 

exp{=R— ^tr [d 3 xe lluX a fl (F uX + A u a x )}. (27) 

107T J 

We can rewrite I [A] in the form 

exp(-/[A]) = J Pa^expj-tr J d 3 x (^S^a u ± ^-*F^)}, (28) 
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with 

S^ ab = ^lM ab T ^ W QC ^a (29) 

and ^ 

= -^ixv\F v x- (30) 

Then, if we perform the path-integral over a M we find 

nA] = T ^- 2 trJd*x*F,S^F u . (31) 
Hence, we see that Zj can be written in the form 

Zj = J VA^MSfcs[A}) = Z FCS (32) 

where 

S FOS [A] = ±-^S cs [A] + tr J d\ *F ll S flv *F v (33) 

Then, using eq. ( |25"1) we have the relation 

Z Th « Z FC5 (34) 

which shows the equivalence (to order 1/m) of the fermionic Thirring model 
and a bosonic theory with action Sfcs- Let us notice that only in the limit 
g 2 — > the action Sires reduces to the Yang-Mills-Chern-Simons action. 

3 Bosonization Identities 

In order to infer the bosonization identities which derive from the equivalence 
found in the last section, we add a source for the Thirring current: 

L T hh\ = L Th + tr J d\fb^ (35) 
Then, instead of (|5|), the partition function now reads 
Z Th [b}= Jv^V4>Va^exp{-Jd 3 x$(i$ + m + 4 + b')^ + 1 ^tra< 1 a^}. (36) 
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Or, after shifting a M — > a M — 6^, 



Z Th [b] = exp{-^tr J d 3 xb^} x J X?a M exp{-S e// [a] + ^ tr / 

^ 9 (37) 

with 5 e //[a] still given by eq. dl2|) . Then, we can again establish to order 1/m 
the connection between the Thirring and self-dual models, in the presence of 
sources: 

Z Th [b] = exp{-^tr J d 3 xb^} x J Va^ exp{-5 5D [a] + ^ tr / d 3 xb^}. 

(38) 

Z T/l [6] = exp{-^ tr J d 3 xb^ 1 } x Z SD [6], (39) 



or 



where 



Z SD [b] = J Va^ exp{-S SD [a] + -j tr | d 3 x & M a"} (40) 

In order to connect this with the vector boson system, let us again consider 
the interpolating Lagrangian Lj (eq. [15]), but now in the presence of sources: 

L I [a,A,b]=L I [a,A] + -£ 1 tia IJ ,U t . (41) 

By shifting as before A^ — > A^ — a M and integrating out A^, one easily shows 
that the corresponding partition function Zj[b] coincides (up to a normaliza- 
tion factor) with Z$D[b], 

Z SD [b] = Zi[b\ (42) 

and then 

Z Th [b] = exp{-^tr jd 3 xb^} x Zj[b}. (43) 
If we integrate in the inverse order we have (after shifting a M — > 2a fJi ), 

Zj[b] = ( VA^expiT^-ScslA}) x exp(-J[A,6]) (44) 

with 

//" I 2 87T 

(45) 
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This means that Zi[b] can be written in the form 
Zj[b] = fvA.expiT-^-ScslADx 

J 107T 

so we have the relation 

Z Th [b] = exp{-^tr J d'xb^x J VA^exp(-S FC s[A])x (47) 

exp{-^-tr fd s xb^S^b v } x eX p{±-^-tr [d 3 xb^S^*F v }. 
8g 4 j lbirg z J 

This is, in its most general form, the result we were after. It provides a 
complete low-energy bosonization prescription, valid for any g 2 , of the matrix 
elements of the fermionic current. Since after differentiating we must set 
6 M = 0, we see that, as suggested by eq. fl34"|), the bosonized version of the 
Thirring model is described by the action Spcs m ec i- (23)- However, the 
bosonization rule for the fermionic current is not simple. For instance, from 
eq. (|47| ) we get, up to contact terms, 



i 



- ± y^ s ^k, (48) 

However, higher derivatives respect to the sources lead to more involved 
bosonic equivalents of the vacuum expectation value of products of fermionic 
currents. 

This more complex structure should not come as a surprise: even in two- 
dimensions a simple bosonization procedure applies only to free fermions [|]]. 
Then, with that in mind, we restrict the discussion that follows to the g 2 ^ 
limit, where with eq. fl29|), eqs. (|33| , [47]) reduce to 

Z Th [b] = JvA ll ex P (T^S cs [A})xex P {±^-tTjd 3 xb^}x 

exp{±^-tr Jd 3 xe^b^[A a ,b u ]}, (49) 

so, in that limit, the ground state fermionic current maps to 



n - (50) 
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Concerning the factor in (pE9| ) which is quadratic in the sources, its first con- 
tribution will arise when computing current- current correlation functions. 
One can see however that the resulting commutator algebra, obtained via 
the Bjorken- Johnson- Low method from these correlation functions, is not 
modified by the quadratic term. In this sense, eq.(^) gives the bosonization 
mapping for non-abelian free fermions in 3 dimensions as the natural gener- 
alization of the result obtained in ref.|l| for the abelian case. Of course the 
fact that our results only hold at long distances makes the analysis of the 
commutator algebra, which tests short distances, not completely reliable. 

We thus see that the non-abelian bosonization of free SU(N) massive 
fermions in 2 + 1 dimensions leads to the (level k = 1) SU (N) Chern-Simons 
theory, with the fermionic current being mapped to the dual of the gauge 
field strength. As stated earlier, this result holds only for length scales large 
compared with the Compton wavelength of the fermion, since our results 
were obtained for large fermion mass. It is important to notice that the limit 
g 2 —>■ to which we restrict henceforth, corresponds to free fermions but not 
to an abelian gauge theory. On the contrary. is the full non-abelian 
field strength (cf. eq. (|16D ) and the Yang-Mills coupling is proportional to 
1/g 2 — > oo, which is why we are left with a pure Chern-Simons theory and 
not with a mixed Yang-Mills- Chern-Simons action. 

Eq. ( [50D gives a natural non-abelian extension of the abelian bosonization 
rule obtained in ref. In the abelian case one can interpret the (2 + 1)- 
dimensional bosonization formula (which is identical to eq. ([5(]) but with 
the abelian curvature) as the analog of the (1 + l)-dimensional result 
ipry^ip _=> (iy \fn)e^ v d v '(j). (The factor of i in the expression for the current in 
eq. (|50D appears because we are working in Euclidean space). To establish 
this correspondence also in the non-abelian model, one should remember that 
in this last case the two dimensional bosonization identity reads [|J 



with fermions in the fundamental representation of the group G, and h an 
element of G. 

To make contact with our (2 + l)-dimensional result, let us first note that 
in the g 2 — > limit the resulting bosonic action (the Chern-Simons action) is 



3+ 



3 




(52) 



(51) 
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gauge invariant and so a gauge fixing is required. The natural gauge in order 

to compare the results in 1 + 1 and 2 + 1 dimensions is the A 3 = gauge. 
Moreover, one can write 

A± = A X ± iA 2 (53) 

in the form 

A+ = -ihr x d+h (54) 

A_ = -if^d.f (55) 

With this, one has from (|50|) , 

J+ - ±^d 3 (h~ 1 d + h) (56) 

3 - T^0 3 (f-^f) (57) 

These are the (2 + l)-dimensional analogs of the two-dimensional formulas 
fl5lD- ([52"|) . Concerning the additional component j 3 , one has 



j 3 = ±JL (d+if'^-f) -d-{h- l d + h) + i{h- l d + hj- l d^f)) (58) 
4 Wilson Loops 

Perhaps the most interesting aspect of the bosonization identities of eqs. (f49l)- 
( |50|) is the promotion of the global 5^7 (AT) symmetry of the free fermions to 
a local gauge symmetry in the bosonic theory. To explore the contents of 
this bosonization rule we consider the natural objects in the gauge theory, 
namely the vacuum expectation value of Wilson loops. In the Chern-Simons 
theory they measure topological invariants determined by the linkings of the 



loops and by the topology of the base manifold |15 |. For one loop V 



W[T] = ti P exp(i £dx» A^) (59) 

where P denotes the path ordering of the exponential, and the trace is taken 
in the representation carried by the loop. According to the bosonization 
prescription of eq. (|50D, to relate this operator to the fermionic theory we 
must express W[T] in terms of the field strength F^ u rather than the potential 
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A,j,. In the abelian case this can always be done by means of Stokes theorem. 
As discussed in ref. 0, this leads to an explicit mapping between abelian 
Wilson loops and non-local fermionic operators. Hence, in this way, the 
latter are related to the linking of loops and thus probe the generalized 
statistics of the external particles that propagate along those loops. One 
way to extend that analysis to the non-abelian case is to use the non-abelian 



extension of Stokes theorem developed in [16|]. For an arbitrary loop T = <9E, 
the boundary of a surface E, one has 

W[dE]=tvP t exp{i dt ds- --H/- 1 [ s S(t) ]F^(S(t, S ))H/[ s E(t)o]} 

Jo Jo as at 

(60) 

Here E is looked upon as a sheet, that is, a one parameter family of paths 
parametrized by t, < t < 1. For each t, E(t) is a path, itself parametrized 
by s, < s < 1, with fixed end-points: <9E(t, s) / 'dt = at s = 0, 1. For a 
given t, s S(t)o denotes the segment of the path E(t) connecting the points 
E(t, 0) and E(t, s), and W[ s £(i)o] is the corresponding (open) Wilson line. 
Finally, Pt in eq. ( |60D denotes ordering of the t integral, while the s integral 
is not ordered (although there is an s-ordering inside each H / [ s E(t) ].) 

In the abelian case the two open Wilson lines W[ s £(£) ] in eq. ( loTf ) cancel 
each other and one recovers the usual Stokes theorem, involving only the 
gauge field strength. In the non-abelian case, however, the factors H^[ s E(t) ] 
are needed for gauge invariance, and introduce an explicit dependence of the 
Wilson loop operator on the gauge potential A^. Thus, as opposed to the 
abelian case, the non-abelian Wilson loop operator cannot be mapped in a 
straightforward way to a fermionic operator through the bosonization rule in 
eq. ®. 

For planar loops this difficulty is only apparent. Indeed, consider W[<9E], 
with E contained, say, in the (1,2) plane. Imposing the A 3 = gauge 
condition, there is a remnant gauge freedom for the Ai and A 2 components 
in the (1,2) plane, which is the symmetry of a 2-dimensional gauge theory in 



that plane. As discussed in |lq] , one can use that gauge symmetry, together 
with the freedom of parametrization of the surface E, so the open Wilson 
line elements in the right hand side of eq. (|K]) become the identity. More 



precisely, choosing the gauge condition A 2 = on the E-plane, eq. (pOf) can 
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be simplified to 

r 1 r 1 dY7 

W[dE]=txP t exp{i dt ds- — F^{E{t,s))} (61) 

Jo Jo OS Ot 

provided that S is parametrized so as to have dY^/dt and dH/ds parallel 
to the xi and X2 axis, respectively. This apparent breaking of rotational 
invariance, which includes the presence of t-ordering but not of s-ordering, 
is a consequence of the A 2 = gauge condition on the S-plane, and will be 
removed by the functional integral over the gauge fields. Then, writing 

J[E] = trP t exp{±47r J dt J ds Qg ^ e^ A j A [S(t, s)}} (62) 

with jp the fermionic current in eq. @, the bosonization formula (|50D gives 

(J[n) ferm =(W[dE])cs (63) 

where in the left hand side the subindex 'ferm' stands for free fermions. 
This is the non-abelian generalization of the result obtained in [0. It relates 
a suitably defined non-abelian flux of the fermionic current through a flat 
surface, and the Wilson loop associated to the boundary of that surface, 
with both quantities in the same representation of the group. 

It should be stressed that the bosonic side of this relation is, by definition, 
independent of the surface S and its parametrization. In the fermionic side, 
however, this is not obvious. The relation was derived assuming a flat surface 
S, and it is tempting to assume that this may be extended to smooth defor- 
mations away from the plane. But more relevant is the apparent breaking 
of rotational invariance in the fermionic side due to the remaining t-ordering 
in eq. (|62|). This should certainly be expected to be taken care of by the 
particular parametrization assumed above for S. Indeed, one should expect 
that the very need of a parametrization and of a matching ordering of the 
surface integral of the fermionic current, is just a limitation of our present 
analysis. In addition, as is well known, the expectation value of the Wilson 
loop is singular and must be regularized. A natural and consistent regular- 
ization scheme is provided by the framing of the loop [15]. In the case of a 



non-intersecting loop on a plane, considered here, that framing can be chosen 
also as a plane loop, not intersecting itself nor the original loop. Again, it is 
not clear at this point how these singularities in the bosonic side will show 
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up in the (free) fermionic side, and what role will the framing play from the 
fermionic point of view. 

It is natural to ask whether this analysis can be extended to several loops 
and their possible linkings, as done in for the Abelian case. In the bosonic 
side one is interested in the expectation value (W / [r 1 ]W / [r 2 ])c5 or, better yet, 
the ratio 

For non intersecting loops this is a well defined, non singular object in the 
Chern-Simons theory, which depends only on the linking of the two loops 
r"i and T 2 Assuming this to be non-trivial (and non-singular), the two 
loops cannot be flat and lying on the same plane, so the previous construction 
fails. But once the ratio (|64"|) has been computed in the Chern-Simons theory, 
we can take the limit in which the two loops collapse onto a single plane. 
This is a singular limit in which the loops necessarily intersect each other. 



Their linking is not well defined any more, and the value of (64) depends 
on the initial non-singular loops used in the computation. However, at the 
classical level, before the functional integral is performed, we can repeat the 
previous construction with no difficulties for any arrangement of loops on the 

(65) 



plane |16 |. Thus, formally we can write 

(J\E X U S 2 ]) /erm {W[dZ x ]W[dY, % ])cs 



{J[Z x \) ferm {J[Z 2 \) }erm {W[dZ x \)cs{W[dZ 2 ]) cs 

where both surfaces Si and S 2 are contained in the same plane. As we just 
stated, the bosonic side of this relation will be ill defined in general. But it can 
be given a well defined meaning by lifting the loops <9£j from the plane to non 
intersecting three-dimensional loops IV This can be done in different ways, 
specifying different linkings of the loops Tj compatible with the intersections 
of their projections <9£j onto the plane. Correspondingly, in the fermionic 
side, the surface Si U S 2 must be complemented with a prescription stating 
the way in which the two surfaces Ej overlap. The different possible liftings 
of the loops specify different overlaps of the surfaces, as ilustrated in Fig. (1). 
In this way, relation fl6"5|) (and its generalizations) can be viewed as a defining 
relation, through bosonization, of the vaccuum expectation value of the flux 
of the fermionic current through surfaces with foldings. 
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5 Conclusions 



We close with a few remarks on the nature of the mapping discussed here. 
In this paper we showed that the low energy sector of the massive SU(N) 
Thirring model in 2 + 1 dimensions is equivalent to the long distance regime of 
a non-abelian gauge theory, closely related to the Yang-Mills-Chern-Simons 
gauge theory. In the weak coupling limit, the two theories become identi- 
cal. It is worthwile to stress that this mapping only holds at long distances. 
In that regime, the gauge theory is a topological field theory and so is the 
fermion theory. Secondly, just as in the abelian case, we discover the ex- 
istence of operators of the Fermi theory which ought to exhibit fractional 
statistics. However, unlike the abelain theory, these objects are substantially 
more complex. Finally, the bosonic theory is, essentially, a level k — 1 SU (N) 
Chern-Simons theory. It would be interesting to find a fermionic analog of a 
Chern-Simnos theory with level higher than one. 
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Figure 1: Different overlaps of the surfaces Si and S2 on a plane, determined 
by the possible liftings of the loops <9£i and <9E 2 away from the plane. 
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